A Generic Pharmacokinetic Model for Quantifying Mother-to-Offspring Transfer of Lipophilic Persistent Environmental Chemicals.

Lipophilic persistent environmental chemicals (LPECs) can accumulate in a woman's body and transfer to her developing child across the placenta and via breast milk. To assess health risks associated with developmental exposures to LPECs, we developed a pharmacokinetic (PK) model that quantifies mother-to-offspring transfer of LPECs during pregnancy and lactation and facilitates internal dosimetry calculations for offspring. We parameterized the model for mice, rats, and humans using time-varying functions for body mass and milk consumption rates. The only required substance-specific parameter is the elimination half-life of the LPEC in the animal species of interest. We used the model to estimate whole-body concentrations in mothers and offspring following maternal exposures to hexachlorobenzene (HCB) and 2,2',4,4',5,5'-hexachlorobiphenyl (PCB 153) and compared these with measured concentrations from animal studies. We also compared estimated concentrations for humans to those generated using a previously published human LPEC PK model. Finally, we compared human equivalent doses (HEDs) calculated using our model and an allometric scaling method. Estimated and observed whole-body concentrations of HCB and PCB 153 in offspring followed similar trends and differed by less than 60%. Simulations of human exposure yielded concentration estimates comparable to those generated using the previously published model, with concentrations in offspring differing by less than 12%. HEDs calculated using our PK model were about 2 orders of magnitude lower than those generated using allometric scaling. Our PK model can be used to calculate internal dose metrics for offspring and corresponding HEDs and thus informs assessment of developmental toxicity risks associated with LPECs.

Framework for risk assessment of PFAS utilizing experimental studies and in-silico models.

Perfluoroalkyl substances (PFAS), especially PFOS and PFOA, are two widely used synthetic chemicals that can impact human health based on evidence from animal and epidemiologic studies. In this paper, we have reviewed and summarized the influence of PFAS exposure on health, pointing the quality of evidence, and applied translational techniques to integrate evidence for PFAS policy making. This is the first review where highly referred articles on PFAS used for policymaking by several regulatory agencies were collected and evaluated based on the review guidelines developed by the US National Toxicology Program's Office of Health Assessment and Translation (OHAT) review guidelines. Several limitations were observed, including co-exposure to multiple chemicals and limited measurement of primary and secondary outcomes related to specific toxicity. However, data from all the studies provided a moderate to strong level of confidence for link between PFAS exposure and different adverse outcomes. Secondly, for translating the risk to humans, an in-silico model and scaling approach was utilized. Physiologically based pharmacokinetic model (PBPK) was used to calculate the human equivalent dose (HED) from two widely accepted studies and compared with tolerable daily intakes (TDIs) established by various regulatory agencies. Inter-species dose extrapolation was done to compare with human the relevance of dosing scenarios used in animals. Overall, a framework for translation of risk was proposed based on the conclusions of this review with the goal of improving policymaking. The current paper can improve the methodological protocols for PFAS experimental studies and encourage the utilization of in-silico models for translating risk.

Tissue distribution study of perfluorooctanoic acid in exposed zebrafish using MALDI mass spectrometry imaging.

Perfluorooctanoic acid (PFOA) as an emerging environmental contaminant, has become ubiquitous in the environment. It is of significance to study bioconcentration and tissue distribution of aquatic organisms for predicting the persistence of PFOA and its adverse effects on the environment and human body. However, the distribution of PFOA in different tissues is a complex physiological process affected by many factors. It is difficult to be accurately described by a simple kinetic model. In present study, a new strategy was introduced to research the PFOA distribution in tissues and estimate the exposure stages. Zebrafish were continuously exposed to 25 mg/L PFOA for 30 days to simulate environmental process. Matrix-assisted laser desorption ionization mass spectrometry imaging (MALDI-MSI) method was used to monitor the spatio-temporal distribution of PFOA in zebrafish tissues. By analyzing the law of change obtained from the high spatial resolution MSI data, two different enrichment trends in ten tissues were summarized by performing curve fitting. Analyzing the ratio of two types of curves, a new "exposure curve" was defined to evaluate the exposure stages. With this model, three levels (mild, moderate, and deep pollution stage) of PFOA pollution in zebrafish can be simply evaluated.

The oral bioavailability of di-2-ethylhexyl phthalate (DEHP), di-isononyl phthalate (DiNP) and di-(isononyl)-cyclohexane-1,2-dicarboxylate (DINCH®) in house dust.

Phthalates and other plasticizers are detected in high amounts in the indoor environment and therefore house dust can be an exposure source. Especially children have a relatively high unintended uptake of house dust, thus a higher exposure to plasticizers compared to adults may be possible. As accurate as possible exposure assessment data of the oral bioavailability of these compounds are necessary, however only one in vivo study with piglets is available so far. The aim of this study was to examine the oral bioavailability of phthalates and DINCH® in humans, which occur in typical house dust samples. We focused on the high molecular weight phthalates DEHP and DINP and their substitute DINCH®. Eleven volunteers ingested 6 g of house dust sieved to 2 mm. The urine was collected over a period of 36 h. The excreted plasticizers metabolites were quantified by an LC-MS/MS method. The mean recovery of urine metabolites was 51 % ± 20 % for DEHP, 26 % ± 13 % for DINP and 19 % ± 6% for DINCH® based on the parent compounds administered as dust samples. The metabolites of DEHP, DINP and DINCH® reached their maximum concentration after 2-19 hours post dose in urine. The bioavailability of DEHP was in agreement among the different dust samples. For DEHP, we were able to confirm previous findings from the oral bioavailability study with piglets and we could not observe a significant difference between the dust particle size (65 µm vs 2 mm) and the bioavailability. Considering the observed bioavailability, an estimated dust intake of 50 mg/d for toddlers can substantially contribute to the total plasticizer exposure.

Transformation of Tetracycline by Manganese Peroxidase from Phanerochaete chrysosporium.

The negative impacts on the ecosystem of antibiotic residues in the environment have become a global concern. However, little is known about the transformation mechanism of antibiotics by manganese peroxidase (MnP) from microorganisms. This work investigated the transformation characteristics, the antibacterial activity of byproducts, and the degradation mechanism of tetracycline (TC) by purified MnP from Phanerochaete chrysosporium. The results show that nitrogen-limited and high level of Mn2+ medium could obtain favorable MnP activity and inhibit the expression of lignin peroxidase by Phanerochaete chrysosporium. The purified MnP could transform 80% tetracycline in 3 h, and the threshold of reaction activator (H2O2) was about 0.045 mmol L-1. After the 3rd cyclic run, the transformation rate was almost identical at the low initial concentration of TC (77.05-88.47%), while it decreased when the initial concentration was higher (49.36-60.00%). The antimicrobial potency of the TC transformation products by MnP decreased throughout reaction time. We identified seven possible degradation products and then proposed a potential TC transformation pathway, which included demethylation, oxidation of the dimethyl amino, decarbonylation, hydroxylation, and oxidative dehydrogenation. These findings provide a novel comprehension of the role of MnP on the fate of antibiotics in nature and may develop a potential technology for tetracycline removal.

Discrimination of active and inactive substances in cytotoxicity based on Tox21 10K compound library: Structure alert and mode of action.

In vitro cytotoxicity assay is an ideal alternative method for the in vivo toxicity in the risk assessment of pollutants in environment. However, modes of action (MOAs) of cytotoxicity have not been investigated for a wide range of compounds. In this paper, binomial and recursive partitioning analysis were carried out between the cytotoxicity and molecular descriptors for 8981 compounds. The results showed that cytotoxicity is strongly related to the chemical hydrophobicity and excess molar refraction, indicating the bio-uptake and chemical-receptor interaction through π and n electron pair play important roles in the cytotoxicity. The decision tree derived from recursive partitioning analysis revealed that the studied compounds could be divided into 25 groups and their structural characteristics could be used as structure alert to identify active and inactive compounds in cytotoxicity. The descriptors used in the decision tree revealed that chemical ionization and bioavailability could affect the cytotoxicity for ionizable and highly hydrophobic compounds. Comparison of MOAs based on Verhaar's classification scheme showed that many inert or less inert compounds were inactive substance, and many reactive or specifically-acting compounds were active substances in the cytotoxicity. In vitro toxicity assay instead of in vivo toxicity assay can be used in the environmental hazard and risk assessment of organic pollutants. The descriptors used in the binomial equation and decision tree reveal that chemical hydrophobicity, ionization and solubility play very important roles for identification of active and inactive compounds. The results obtained in this paper are valuable for understanding the modes of action in cytotoxicity and in vivo-in vitro toxicity relationship.

The Role of Biotransformation in the Activity of Endocrine Disruptors.

An "endocrine disruptor" has been broadly defined as an exogenous chemical that interferes with the production, release, transportation, metabolism, binding, action, or the elimination of endogenous hormones, which are responsible for homeostasis, reproduction, development or behaviour. Diverse groups of chemicals such as pharmaceuticals, phytoestrogens, natural hormones, and synthetic chemicals such as pesticides, plasticizers, phthalates, parabens, polychlorinated/polybrominated biphenyls, bisphenols are shown to interfere with the endocrine system, and they have been defined as EDs in the last three decades. As for all chemicals, the biotransformation of EDs has a decisive role in their potential toxic effects. Humans are exposed to vast amounts of diverse chemicals throughout their lives. Fortunately, most of the chemicals are converted via biotransformation reactions catalyzed by the enzymes, into more hydrophilic metabolites, which are readily excreted in urine or bile. Biotransformation reactions resulting in less toxic metabolites are known as detoxification. However, some biotransformation reactions are called bioactivation, in which more toxic metabolites are formed. In the case of EDs, metabolites formed via bioactivation usually have a higher affinity for a hormone receptor or induce/inhibit an enzyme involved in the synthesis or catabolism of an endogenous hormone more dramatically compared to their parent compound. In the present review, the role of bioactivation in endocrine modulating effects of chemicals from all groups of EDs, namely endogenous estrogens, phytoestrogens, synthetic/industrial chemicals, and pharmaceuticals it can be were discussed.

Effects of curcumin on the bioavailability of dioxin-like pollutants in rats.

The effects of curcumin on the bioavailability of polychlorinated dibenzo-p-dioxins/furans (PCDD/Fs) and dioxin-like polychlorinated biphenyls (DL-PCBs) were investigated in Sprague-Dawley rats. Tetra- and penta-chlorinated PCDFs had the lowest bioavailability and hexa-chlorinated PCDD/Fs had the highest, while there was no obvious change in that of DL-PCBs. Curcumin markedly reduced the toxic equivalent (TEQ) of PCDD/Fs in rats, illustrating the potential to competitively inhibit absorption of PCDD/Fs by the epithelial cells of the small intestine due to the similar chemical structure (diphenyl) between curcumin and PCDD/Fs. Moreover, curcumin lowered the TEQ of DL-PCBs in the liver of male rats, but not female rats. The significant decrease in the bioavailability of PCDD/Fs and DL-PCBs demonstrates the potential detoxification mechanisms of curcumin.

THE DEVELOPMENT OF A HIGH-RESOLUTION MASS SPECTROMETRY METHOD FOR ULTRA-TRACE ANALYSIS OF CHLORINATED DIOXINS IN ENVIRONMENTAL AND BIOLOGICAL SAMPLES INCLUDING VIET NAM ERA VETERANS.

Chlorinated dioxins are labeled and recognized by both the World Health Organization and the United Nations Environmental Programme (UNEP) as "persistent organic pollutants". Their potential for high toxicity is one of the primary factors behind intense public and regulatory scrutiny and the need to measure the compounds at very low limits, specifically the isomer 2,3,7,8-tetrachlorodibenzo-p-dioxin (2,3,7,8-TCDD). This article highlights the early mass spectrometry methods to investigate, detect, confirm, and quantify chlorinated dioxins and the initial applications involving human biomonitoring, as attempts were made to attribute health effects to TCDD exposure. This effort represented a complex and difficult scientific response to the pressing need to investigate expected exposures and alleged subsequent medical effects, which in the case of the Viet Nam veterans was being attempted a decade or more after their exposure. It is noteworthy that this method and its development touched on delicate issues involving human subjects, war veterans, environmental contamination, and was difficult not only scientifically, but for ethical and political reasons as well. Stable-isotope dilution with analysis by gas chromatography/high-resolution mass spectrometry (GC/HRMS) became the method of choice because of its ability to monitor characteristic ions and isotope ratios to quantify and qualify/confirm the analyte in the presence of coextracting and coeluting interferences at these low levels with the highest possible confidence. This method was rigorously tested and validated before it was used to discover and monitor levels in the environment and in various populations at then unprecedented low levels. These early studies demonstrated the feasibility of monitoring dioxins in humans even decades after exposure, and led to the detection of 2,3,7,8-TCDD in the general population as well as specific overexposed populations. These studies also provided strong evidence regarding the origins of the 2,3,7,8-isomer in the environment. © 2020 John Wiley & Sons Ltd. Mass Spec Rev.

The Role of Metabolism in Toxicity of Polycyclic Aromatic Hydrocarbons and their Non-genotoxic Modes of Action.

Polycyclic aromatic hydrocarbons (PAHs) represent a class of widely distributed environmental pollutants that have been primarily studied as genotoxic compounds. Their mutagenicity/genotoxicity largely depends on their oxidative metabolism leading to the production of dihydrodiol epoxide metabolites, as well as additional metabolites contributing to oxidative DNA damage, such as PAH quinones. However, both parental PAHs and their metabolites, including PAH quinones or hydroxylated PAHs, have been shown to produce various types of non-genotoxic effects. These include e.g., activation of the aryl hydrocarbon receptor and/or additional nuclear receptors, activation of membrane receptors, including tyrosine kinases and G-protein coupled receptors, or activation of intracellular signaling pathways, such as mitogen-activated protein kinases, Akt kinase and Ca2+-dependent signaling. These pathways may, together with the cellular DNA damage responses, modulate cell proliferation, cell survival or cell-to-cell communication, thus contributing to the known carcinogenic effects of PAHs. In the present review, we summarize some of the known non-genotoxic effects of PAHs, focusing primarily on those that have also been shown to be modulated by PAH metabolites. Despite the limitations of the available data, it seems evident that more attention should be paid to the discrimination between the potential non-genotoxic effects of parental PAHs and those of their metabolites. This may provide further insight into the mechanisms of toxicity of this large and diverse group of environmental pollutants.

Metal remediation potential of naturally occurring plants growing on barren fly ash dumps.

The present study aimed to elucidate the remediation potential of visibly dominant, naturally growing plants obtained from an early colonized fly ash dump near a coal-based thermal power station. The vegetation comprised of grasses like Saccharum spontaneum L., Cynodon dactylon (L.) Pers., herbs such as Tephrosia purpurea (L.) Pers., Sida rhombifolia L., Dysphania ambrosioides (L.) Mosyakin & Clemants, Chromolaena odorata (L.) King & H.E. Robins along with tree saplings Butea monosperma (Lam.) Taub. The growth of the vegetation improved the N and P content of the ash. Average metal concentrations (mg kg-1) in the ash samples and plants were in order Mn (345.1) > Zn (63.7) > Ni (29.3) > Cu (16.8) > Cr (9.9) > Pb (1.7) > Cd (0.41) and Cr (58.58) > Zn (52.74) > Mn (39.09) > Cu (10.71) > Ni (7.45) > Pb (5.52) > Cd (0.14), respectively. The plants showed fly ash dump phytostabilization potential and accumulated Cr (80.19-178.11 mg kg-1) above maximum allowable concentrations for plant tissues. Positive correlations were also obtained for metal concentration in plant roots versus fly ash. Saccharum spontaneum showed highest biomass and is the most efficient plant which can be used for the restoration of ash dumps.

Effects of inhaled combined Benzene, Toluene, Ethylbenzene, and Xylenes (BTEX): Toward an environmental exposure model.

Combined environmental exposures to the volatile organic compounds (VOCs) Benzene, Toluene, Ethylbenzene, and Xylene (BTEX) pose clear risks to public health. Research into these risks is under-studied even as BTEX levels in the atmosphere are predicted to rise. This review focuses on the available literature using single- and combined-BTEX component inhaled solvent exposures in animal models, necessarily also drawing on findings from models of inhalant abuse and occupational exposures. Health effects of these exposures are discussed for multiple organ systems, but with particular attention on neurobehavioral outcomes such as locomotor activity, impulsivity, learning, and psychopharmacological responses. It is clear that animal models have significant differences in the concentrations, durations and patterns of exposure. Experimental evidence of the deleterious health and neurobehavioral consequences of exposures to the individual components of BTEX were found, but these effects were typically assessed using concentrations and exposure patterns not characteristic of environmental exposure. Future studies with animal models designed appropriately to explore combined BTEX will be necessary and advantageous to discovering health outcomes and more subtle neurobehavioral impacts of long-term environmental exposures.

Accumulation of trace metals in eggs and hatchlings of Chelonia mydas.

BACKGROUND: The objective of this study was to verify the accumulation of trace metals in eggs and hatchlings of Chelonia mydas, evaluating if metal accumulation is originated from maternal transfer and/or from the incubation environment. Other assessments were also performed, as metal distribution in different tissues (blood, kidney, liver, muscle, and turtle shells) of newly hatched turtles, and genotoxic analysis, to verify possible damages caused by the presence of metals. METHODS: The assessments were carried out by quantifying Cd, Ni, Pb, Mn and Fe in egg sample collected during laying time (eggshells (ELT) and egg content (EC)), eggshells from newly hatched turtles (ENH), hatchlings tissues (H - blood, kidney, liver, muscle, and shell)) (n = 18 for each biological sample - 3 of each nest) and nest sediments (n = 6, one of each nest). Comparative analysis were made between ELT and ENH, as well as between egg content (EC) and the sum of tissue samples from hatchlings, using Mann-Whitney hypothesis test (p < 0,05). The amount of metals in different hatchling was quantified and followed by the Dunn post-test. A principal component analysis (PCA) was also employed. RESULTS: Metals studied were found in all investigated samples. The concentration of a great amount of investigated metals was significantly higher (P=<0.001) in eggshells from ENH than in ELT. An increase in Cd (2.16-fold), Pb (3.47-fold), Fe (6.83-fold) and Mn (195.57-fold) concentration was noticed in ENH. We also observed an increase in Fe (1.59-fold), Mn (1.74-fold) and Ni (1.59-fold) concentration in hatchling, when compared with EC, due to transfer from nest sediments. In relation to the hatchling's tissues, blood was shown to accumulate higher concentrations of Ni and Pb, while shells accumulated more Cd and Fe, and Mn is more associated with liver and kidney. Fe was the highest accumulated metal in both tissues, and muscles presented discrete concentrations of Ni, Mn, and Pb. A mean concentration of 1.25‰ MN was obtained in C. mydas hatchlings, indicating that the accumulation of metals in hatchlings didn't cause toxicology effects. CONCLUSION: Hatchlings accumulate metals through the maternal and sediment transfer, although the levels of metal accumulation were not enough to cause genotoxic damage.

Body Representations of Internal Pollution: The Risk Perception of the Circulation of Environmental Contaminants in Pregnant and Breastfeeding Women in Spain.

In this article, we analyze how pregnant and breastfeeding women perceive the inside of their bodies as well as their thoughts regarding the accumulation and elimination of chemical compounds present in food, and how these are then transmitted to the fetus. We explore different social perceptions of risk regarding the circulation of chemical compounds inside the body using qualitative research based on the technique of body mapping, comprised of women's drawings of their bodies in combination with comments on the drawings, food diaries and narratives from in-depth interviews. We examine how these 41 women (21 pregnant and 20 breastfeeding) perceive the body's internal mechanisms during the stages of pregnancy and breastfeeding, as well as the circulation of chemical contaminants within it. The body mapping technique allowed us to analyze participants' knowledge of internal pollution, a little-understood process in society. Thanks to these pregnant and breastfeeding women, who made an effort to represent and reflect on these new risks, this study shows that scientists and obstetricians need to collaborate with women in order to better understand and publicize the risks of internal pollution.

Development of a physiologically based pharmacokinetic model of diisononyl phthalate (DiNP) in pregnant rat and human.

A physiologically based pharmacokinetic (PBPK) model for di-isononyl phthalate (DiNP) was developed by adapting the existing models for di(2-ethylhexyl) phthalate (DEHP) and di-butylphthalate (DBP). Both pregnant rat and human time-course plasma and urine data were used to address the hydrolysis of DiNP in intestinal tract, plasma, and liver as well as hepatic oxidative metabolism and conjugation of the monoester and primary oxidative metabolites. Data in both rats and humans were available to inform the uptake and disposition of mono-isononyl phthalate (MiNP) as well as the three primary oxidative metabolites including hydroxy (7-OH)-, oxo (7-OXO)-, and carboxy (7-COX)-monoisononyl phthalate in plasma and urine. The DiNP model was reliable over a wide range of exposure levels in the pregnant rat as well as the two low exposure levels in humans including capturing the nonlinear behavior in the pregnant rat after repeated 750 mg/kg/day dosing. The presented DiNP PBPK model in pregnant rat and human, based upon an extensive kinetic dataset in both species, may provide a basis for assessing human equivalent exposures based upon either rodent or in vitro points of departure.

Autophagy role in environmental pollutants exposure.

During the last decades, the potential harmfulness derived from the exposure to environmental pollutants has been largely demonstrated, with associated damages ranging from geno- and cyto-toxicity to tissue malfunction and alterations in organism physiology. Autophagy is an evolutionarily-conserved cellular mechanism essential for cellular homeostasis, which contributes to protect cells from a wide variety of intracellular and extracellular stressors. Due to its pivotal importance, its correct functioning is directly linked to cell, tissue and organismal fitness. Environmental pollutants, particularly industrial compounds, are able to impact autophagic flux, either by increasing it as a protective response, by blocking it, or by switching its protective role toward a pro-cell death mechanism. Thus, the understanding of the effects of chemicals exposure on autophagy has become highly relevant, offering new potential approaches for risk assessment, protection and preventive measures to counteract the detrimental effects of environmental pollutants on human health.

The influence of five metallic nanoparticles on the expression of major drug-metabolizing enzyme genes with correlation of inflammation in mouse livers.

Metallic nanoparticles (NPs) are widely used in medical preparations. The present study aims to find out the influence of widely used five metallic NPs on the expression of major hepatic drug-metabolizing enzyme (DME) genes. Six groups of BALB/C mice, 7 mice each, were exposed to: Gold NPs, silver NPs, copper oxide NPs, silicon dioxide NPs and zinc oxide NPs, for 21 days. Liver biopsies from all mice were subjected to mouse cyp3a11, cyp2c29, ugt2b1 and interleukin-6 (il6) gene expression quantification using real-time polymerase chain reaction, in addition to inflammatory cell infiltration examination. All tested NPs caused a sharp and significant (ANOVA, p value <0.05) downregulation in the expression of DME genes, with the highest influence was observed in mice exposed to copper oxide NPs. Additionally, all NPs induced hepatic inflammation and upregulated the expression of il6 gene, which were inversely correlated with the expression of DMEs. It is concluded that all tested NPs downregulated the expression of DME genes, with the highest influence exhibited by copper oxide NPs, in correlation with inflammation and il6 gene induction in the liver. Further studies are needed to find out the effect of anti-inflammatory compounds against the alterations induced by metallic NPs exposure on hepatic DMEs.

Ilex Paraguariensis exposition to As and Cd in a closed soilless system.

The exposition of mate (Ilex Paraguariensis A. St.-Hil.) to As and Cd was investigated in plants derived from young mini-cuttings. Mate plants were cultivated in a closed soilless system, composed of coarse sand as substrate and flood fertirrigation. Plantlets were fertirrigated with nutritive solution and As and Cd solutions were added to the nutritive solution in the final concentration of 8 and 17 mg L-1 (As) and of 17 and 33 mg L-1 (Cd) during 14 days. Results show that stem diameter and Dickson quality index (DQI) variables could not be used as a potential indicator of accumulation of As and Cd. The shoot height, number of leaves and chlorophyll index are variables easy and quick to measure and they can be used as parameters to evaluate the stress caused in mate plants cultivation in a closed soilless system. The highest concentration of As and Cd was in roots of plants. Beyond the roots, As and Cd also can be translocated to the leaves achieving high concentrations. In addition, leaves from the treated mate plants were submitted to a hot infusion extraction in order to simulate the traditional beverage and As and Cd were determined in the infusion. Regarding to the infusion procedure, considerable As and Cd amounts were extracted from the leaves leading to conclude that this way of consumption can be an important source of toxic elements for the human diet.

Toxic effects of different-sized graphene oxide particles on zebrafish embryonic development.

Graphene oxide (GO) has broad application potential in many fields, such as biomedicine and energy. Due to the wide-ranging GO applications, its entry into the environment is inevitable along with the potential for ecological and environmental risks. In the present study, we systematically investigated the dose-dependent effects of three different-sized GO particles (50-200 nm, <500 nm, and >500 nm) on zebrafish during the very early developmental stages (4-124 h post-fertilization). The results showed that GOs could accumulate in the eyes, heart, yolk sac, and blood vessels of fish larvae. Consequently, their effects on multiple toxic endpoints were observed, including delayed hatching times, shortened body lengths, alterations in heart rate and blood flow, changes in swimming activity and responses to photoperiod stimulation, and the enhanced activity of total superoxide dismutase, inducible nitric oxide synthase, acetylcholinesterase, caspase-3, and induction of apoptosis-related gene expression. As a result, the occurrence of oxidative stress and the induction of apoptosis are suggested in fish larvae exposed to all three different-sized GO particles. In addition, our results highlight the impacts of waterborne-GO exposure on zebrafish during early development, which were not merely dependent on GO concentration but also on the associated GO sizes. This study hereby provides a basis for the potential ecological and health risks of GO exposure.